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ABSTRACT 



The usefulness of a barotropic prediction model for providing real- 
istic displacements at levels other than 500 mb is suggested to circumvent 
difficulties with current baroclinic models attempting to predict develop- 
ment as well. 

A prediction equation is obtained by making use of concepts from 
Charney*s equivalent barotropic model, ®Arnason B s stratified model, and 
Cressman’s semi-empirical barotropic modelo Further, a stream-function 
prediction scheme, making use of geostrophic and non-dive rgent wind rela- 
tionships, is used to make a series of 24- and 48-hour forecasts of 300-mb 
heights o The vorticity equation used is 



<■ v2 - K 'f fej Tt + J< v, K,P- V) = 0 



where Kj_, K 2 , and are empirical constants 0 Results are presented for 
the effects of adjusting K 2 and 
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1. Introduction. 



Up to the present time, most attempts to predict the pressure height 
at levels other than 500 mb have been made with baroclinic models. The 
results have been disappointing because of inadequately controlled devel- 
opment mechanisms. Observing the persistent skill of barotropic forecasts 
at 500 mb, ‘Arnason [1] suggested the possiblity of devising a barotropic 
prediction model to be applied to other levels, "capable of realistic dis- 
placements of pressure systems. » .but containing no mechanism for develop- 
ment . " 

This investigation introduces an empirically adjusted divergent baro- 
tropic model for application to the 300-mb surface. Concepts from Charney's 
[2] equivalent barotropic model and "Arnason's stratified model [1] are 
incorporated in addition to empirical results by Cressman [33 on control 
of very long atmospheric waves. 
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2. Background. 

Looking first to the equivalent barotropic model, it is recalled that 
partial baroclinicity is permitted in that wind speed, but not direction, 
is allowed to vary with height at any point. If the integrated mean of 
this wind with respect to pressure is determined, the wind, V , and vor- 
ticity, g , at any level can be expressed as a scalar factor. A, of this 
mean. For further simplicity, A is considered to be a function of pressure 
alone . Thus 

V=A( P )V, , (2.D 



where the cap indicates the vertically integrated mean* Climatologically, 
a 

V is observed to occur at about 600 mb e 

Substituting the expressions (2 0 l) in the approximate vorticity 
equation 



H 



+• v- v(£ <•{) 



and integrating with respect to pressure gives 
/\ 



b£ 



4* 



A , A A r \ 

V • V(A a t + ^ " O. 



( 2 . 2 ) 



Here f is the coriolis parameter, and A a is found empirically to have a 

/V X 

value of about 1.25. Multiplying (2.2) by A and substituting A Q- £ and 
A V= V forms the basis for barotropic forecasts at the equivalent baro= 
tropic level . (p*~) s occurring in nature at about 500 mb. 

/\ /\ 

If, however, the expressions (2.1) are substituted for V and £ into 
(2.2), the following is obtained; 



-ff + V + f) =0 



(2.5) 



2 



/\ 



A a 

where For the wind field (2,1), equation (2*3) applies to any 

A 

level with the K 2 appropriate to that level. 

In applying the theory of the stratified model to forecasting, 

1 Arnason arrives at an expression for divergence 



v ' v = -fs(H + > (2,4) 

where J is geopotential, V is the basic flow, and H is the thickness of 
the stratified layer. The parameter, g* , is related to gravity, although 
much smaller, and depends upon static stability and the vertical extent 
of the stratified layer. This expression for divergence is introduced into 
the vorticity equation 

-|^+V-VY) a + V) V • V - O ( 2 -5) 

to give 

+ vv§)=o, (2 - 6) 

where !)=£ + £ is absolute vorticity; T) a is defined for our purposes as 

T) = K-.? < of (2,3). 

The prediction scheme consists of introducing a stream function, y , 
to represent the wind and vorticity 



V = !k x VyJ , X = (2.7) 

where is the unit vector in the vertical. The stream function is ini- 
tially obtained from the balance equation 

y tS ( ) + vy vf - =0 . (2.8) 

Next, in order to simplify the procedure, the quasi-geostrophic relation- 
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ship 



$ = fv <2 - 9) 

is introduced into the divergence tern. The result of substituting (2.7) 
and (2.9) into (2,6) is 

v*y + Kv.'tya) - ^[37 + JCV.V)] =o 

or 

K'h'V) =0, (2,10) 

where y is the stream function corresponding to the basic flow. Equation 
(2.10) is a prognostic equation in a single variable, 9^ . 
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3. Procedure, 



Application of 'Arnason's stratified model involves an interpretation 
of the parameters, y , g', and H. Also, the introduction of a factor, K^, 
into the tendency portion of the divergence term, and a factor, K-j, into 
the advection portion of the divergence tern: permits adjustment of the mag- 
nitude of these terms for various levels. 

The basic flow for this prediction model is taken as a spaced-averaged 
flow at 300 mb, employing a 15-pass, 4-point smoother over the finite-dif- 
ference grid used for computations (grid distance equal to 381 km at 60N 
latitude) : 



yo = 4-(v. •*- u + Vi (3 * 1} 

The subscript zero denotes the central point of a diamond array and the 
subscripts one through four denote the surrounding four grid points , 

In arriving at an interpretation of g’, reference is made to Cressman's 
[3j treatment of a model including a dependence upon static stability. He 
introduced a stability parameter, jx. , which as applied here amounts to a 
substitution of 



yU. „ _l__ 

3 r 



(3.2) 



In a series of 500-mb forecasts Cressman used various values of in his 
divergent barotropic model. It was concluded that yt = 4 affords the de- 
sired control on the ultra-long waves, and gives better forecasts than a 
former Joint Numerical Weather Prediction Unit forecast model. Use has 
been made of this value by the substitution of (3.2) into (2.10). 

In addition, the mean value fbl^has been used for (similar to 

Tim T 

Cressman), taking £500 as the U„ S. standard atmosphere 500-mb height. 
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and f m = f at 45N latitude. 

With these last adjustments, (2.10) becomes 



In finite difference form, (3.3) becomes 



+ I - K 3 .^p^ o 0 ( 9.V)=O , (3.4) 

where SJ7 a , Z4 , and ^ are the finite -difference equivalents of the Lapla- 
cian, differential, and Jacobian operators, respectively; and ^ is the 
sea-level distance between grid points. Forecasts are made in one-hour 
time steps, using a forward difference for the first step and centered 
differences for successive steps. 

Equation (3.4) was programmed for the Control Data Corporation Model 
1604 computer using U, S, Fleet Numerical Weather Facility's 1977-point 
octagonal grid, which covers the Northern Hemisphere from pole to about 
ION. 



Step one in the forecasting procedure consists of solving the balance 
equation for an initial stream function, using analyzed D-values (deviation 
in height from standard). A technique formulated and programmed for the 
CDC 1604 by 'Arnason and Reese at the Fleet Numerical Weather Facility was 
used. It was found that for most of the 300-mb data available for this 
investigation, this method would not converge to the cut-off criterion 
(all residuals less than 0.5 feet in the relaxation procedure as normally 
used at lower levels). In order to proceed with the work, an arbitrary 
cut-off was introduced if non-convergence became apparent. This intro- 
duced errors at a few points of the grid due to residuals which were now 
of the order of three to five feet. These were considered to be insignif- 
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icant with respect to the prediction scheme • 



Step two was solution of equation (3*4) > and this was repeated on 
each day’s data using various values of the factors K 2 and K^* The factor 
i^as taken equal to one and not adjusted in this invest igation 0 Values 
of K 2 and were chosen as work progressed so as to reduce forecast er- 
rors to a minimum. 

Step three consisted of verifying the forecasts against the prepared 
analyses by a simple point -by-point difference, and calculation of the 
pillow and root-mean-square error (RMSE)s 




x 



RMSE 




(3.5) 



x = 1977 points, A = forecast map, B = analysis map. 
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4. Results. 

Table 4.1 lists the results of adjusting the factor K 3 with K 2 held 
equal to one. Numbers to the left of each column are pillow while the 
numbers to the right are RESE. Upper numbers of each paired set are for 
24-hour forecasts and lower numbers are for 48 -hour forecasts. 

Runs were not made for all included values of K3 but minimum errors 
were sought. The numbers in parentheses are values taken from faired 
curves through the values that were run. Fig. 4.1 shows the curves for 
24-hour forecasts; those for 48-hour forecasts would be similar. The 
trends of errors for the successive values of K3 were observed to be quite 
regular, thus justifying an approximation of intermediate values where 
necessary for the averages over all days. In cases where actual reversal 
of error trends was not observed, the fact that differences between suc- 
cessive values of K3 were reduced to zero or one foot was taken as an in- 
dication that reversal was imminent, and adjustment- of K 3 were terminated. 

In one case (17 January) a forecast was made using a positive value 
of K3, showing an error much larger than the optimum for this day. Other 
work on this term not presented here further confirms the large errors for 
values of K3, including K3 equal zero (exclusion of the term), beyond the 
range indicated herein. 

The column entitled "Balancing errors" indicates the differences be- 
tween (a) the height fields obtained after balancing to the stream func- 
tion and inverting again to height, and (b) the initial height field. 

Table 4.2 lists the results of adjusting the value of K 2 , Since K 3 
was first adjusted, optimum values of this factor were used in adjusting 

k 2 . 

It will be noted from table 4.1 that a considerable portion of the 
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forecast errors is apparently due to the stream-function conversion ♦ On 
13 January, the conversion error amounts to Ul% of the total forecast er- 
ror, while on other days the ratios of errors range from 17% to 25 ^. 

The trends of errors for varying all indicate a minimum error in 
the range -0.5 to -1.75 and also indicate the advantage of including 
? Arnason ! s divergence term. From extrapolation of the curves in fig. 4.1, 
errors for equal zero are estimated to average about 20 feet above the 
least average for 24-hour forecasts. 

The use of the factor K 2 does not seem to justify any deviation from 
unity, with improvements ranging from zero to only three feet for 24-hour 
forecasts * 
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Results of adjusting Ko Table 4.1 
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g. 4.1-Trends of RMSS values, 24-hour forecasts. 
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Figs. 4.2 through 4.4 depict the 13 January analyzed height field, 
the 14 January analyzed height field, and the 13 January, 24-hour forecast 
height field verifying on 14 January, respectively. In fig. 4.4 the ini- 
tial (13 January) and observed (14 January) characteristics of certain of 
the systems are superimposed upon the forecast (for 14 January) field to 
indicate contrasting features. Initial features are dotted, and observed 
are dashed. 

The following characteristics may be noted: 

(1) In the western Pacific, the forecast of the short-wave trough 
was very much in agreement with its actual movement and intensity. No 
cyclogenesis occurred in this system and forecast errors were near zero 
in the area of the trough. 

(2) The eastern Pacific low occurring on 13 January almost completely 
diss .pated in 24 hours. The forecast position of this low is nearly that 
of the 14 January trough identifying the system. Height errors in excess 
of 700 feet are due to cyclolysis, which of course is not provided for in 
the model. 

(3) In the alaskan ridge, cyclogenesis occurred from 13 to 14 Janu- 
ary in the NW portion, accounting for errors up to 600 feet. The ridge, 
however, was well forecast. 

(4) The trough off the Atlantic coast of the U. S, was forecast to 
move in agreement with observation, but cyclogenesis again resulted in er- 
rors of about 400 feet. 

(5) The low initially over the central Mediterranean moved into a 
cut-off position south of a blocking high. The forecast movement was in 
fair agreement, although a bit short, but accompanying cyclogenesis raised 
errors to about 600 feet. The forecast position of the blocking high was 
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good, but again, anticyclogenesis in the high and the ridge to the north- 
east introduced forecast errors of 400 feet and 600 feet, respectively. 

Elsewhere on the map, the significant errors seem generally to be due 
to development processes and not to movements of the systems. These ob- 
servations were corroborated by investigation of another sequence of maps 
produced in the investigation. 
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5* Conclusions, 

It is noted from the averaged errors in table 4*1 that there are only 
limited differences in results over the range of values indicated 0 It 
was apparent* however* that sizable differences occur as is extended 
beyond this range* and of particular interest is the advantage over exclu- 
sion of the term (K^ equals zero) . These results indicate an optimum val- 
ue for of about -1*0. 

Even though varying K 2 did not yield significant differences in this 
investigation* it is possible that a greater number of cases might do so, 
particularly if some other season is considered „ Only data from December, 
1961 and January* 1962 were used* 

Taking note of the large errors encountered in converting to a stream 
function, it is possible that a geostrophic model, not requiring such con- 
versions, might reduce total errors somewhat. It is hoped that this in- 
vestigation can be accomplished in the near future. 

The model developed in this paper seems to fulfill its intended pur- 
pose* that is, forecasting realistic displacements unaffected by uncon- 
trolled developments. With an understanding of its limitations, a fore- 
caster should have here a useful tool, when supplemented by techniques for 
development prediction* to forecast pressure systems at the 300-mb level 
more easily and accurately. 
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